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Abstract 

Ideas from causal set theory lead to a fluctuating, time dependent 
cosmological-constant of the right order of magnitude to match cur- 
rently quoted "dark energy" values. Although this effect was predicted 
some time ago [1, 2], it is only more recently that a more detailed phe- 
nomcnological model of a fluctuating A was introduced and simulated 
numerically [3] . In this paper we continue the investigation by study- 
ing the sensitivity of the model to some of the ad hoc choices made in 
setting it up. 

1 Introduction 

As explained in reference [1], a heuristic argument from causal set theory, 
leads one to expect fluctuations in the cosmological constant A which scale 
inversely as the square-root of the spacetime volume. (The reasoning rests on 
the fundamental hypothesis of spatio-temporal discreteness on one hand, and 
on the quantal conjugacy 1 or "uncertainty relation" between four-volume 
and A on the other hand.) Assuming that the fluctuations are centered on 
zero, one obtains for the current epoch the prediction A ~ ±lCP 120 ppz anc fc, 
where ppianck is the Planck density. It is thus natural to interpret the 
observations pointing to an accelerating Hubble expansion as a confirmation 
of this prediction, since the effective energy density corresponding to the 
acceleration is of the same order of magnitude as the predicted fluctuations. 
(The sign of the fluctuations is purely random in this scenario, though with 
a slight bias toward positive values, as described below.) 

Not only does the current value of A receive a natural explanation in 
this way, but computer simulations of a simple phenomenological model [3] 
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for the time-dependence of the fluctuations have confirmed the suggestion 
that a type of "tracking" behavior arises automatically, in the sense that the 
fluctuations would have been comparable to the ambient matter density not 
only now, but in the past as well (and also in the future, for as long as the 
expansion continues). In this way the so called "Why now?" puzzle is also 
resolved. (This puzzle has of course also called forth many other proposed 
solutions, for example [7, 8, 9]. However most of them suffer from the need 
for "fine tuning".) 

The model put forward in [3] derives A (which can be interpreted quite 
generally as the action S of free spacetime per unit 4-volume) from a sum 
of random contributions to S coming from the causal set elements within 
the past light cone of any given spacetime location. It contains a single 
phenomenological parameter which reflects both the magnitude of the in- 
dividual contributions to S, and the conversion factor between spacetime 
volume and number of causal set elements, and which in the absence of 
"fine tuning" would have a value of order unity. Following [3] , we will refer 
to this parameter, which cannot at present be obtained from first principles, 
as a. 

How does this model produce a A having mean zero and fluctuations 
of the desired magnitude? Owing to the random signs of the individual 
contributions to S, their sum will vanish on average, but there will remain 
residual fluctuations in S of order ay/N. Now according to one of the 
fundamental assumptions of causal set theory, the volume of a spacetime 
region must be identified — modulo Poisson fluctuations — with the number 
of elements constituting that region, i.e. V ~ N in natural units. Thus, the 
fluctuations in S produce, at a given cosmic time, a cosmological term with 
a typical magnitude given, as desired, by ay/N /V « a\/V /V « a/W, a 
time-dependent value which diminishes as the volume of the past grows. 2 
If we crudely identify V with (i? -1 ) 4 , H being the Hubble parameter, then 
we can see with the help of the Friedmann equation that one should expect 
the magnitude of the fluctuations in A to track the total energy density: 
A ~ \/l/V ~ 1/ ^(H^ 1 ) 4 ~ H 2 ~ Perineal ■ And exactly this behaviour was 
observed in the numerical simulations of [3]. 

Besides the restriction to spatial homogeneity, the model of [3] contains 
a second ad hoc element of importance. To understand where it comes from, 
recall that the cosmological term A in the Einstein equation, 

-G ab + Ag ab = T ab , (1) 

K 

has classically to be a true constant unless T ab fails to be conserved (or re 
fails to be a constant). 3 In order to interpret the predicted time dependence 

2 The fluctuations are constrained by hand to be spatially homogeneous, this being 
probably the most serious limitation of the model as it has been developed so far. 
3 Here k = 8tvG, where G is Newton's constant. Henceforth we take k — 1. 
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we must therefore depart from the classical field equations. To that end, 
we first assume the universe to be homogeneous and isotropic, so that (1) 
reduces to a set of two ordinary differential equations in a single dependent 
variable, the cosmological scale factor a. For the case of a spatially flat 
universe, which seems to be favored by the cosmological data, these two 
equations are [10] 

3H 2 -A = p (2) 

and 

2a a 2 

— + — -A = -p, 3 

a cr 

where a dot denotes differentiation with respect to proper time r, H = a/a 
is the Hubble parameter, p is the energy density of non-gravitational matter 
(including "dark matter"), and p is its pressure. 

For future reference, notice that if, as is often done, we interpret —g ab A 
as an effective addition to the stress-energy T ab of perfect fluid form, then 
(2) and (3) entail the identifications, p/\_ = A and p^ = —A, corresponding 
to the "fluid equation of state", p\ = — pa. 

Of our two equations for a, the first one or Friedmann equation is only of 
first order in time. It is therefore a constraint equation in the standard ter- 
minology, and we will refer to it in this paper as the Hamiltonian constraint 
equation or HCEq. The second equation is of second differential order, and 
we will refer to it as the acceleration equation or AccEq. 

Now we wish to allow A to vary with time. Since doing so renders (2) 
and (3) incompatible, we can choose at most one of these equations as our 
dynamical guide. In reference [3] we let HCEq play this role, but that is not 
the only possibility, albeit it is a natural choice and it has the nice feature 
that one can interpret the resulting scheme as a local change to the "equation 
of state of A" . 4 Nevertheless, one would like to know how a different choice 
would affect our conclusions. What if we based the model on equation (3) 
rather than equation (2), or on some linear combination of the two, such as 
the trace of the Einstein equation (TrEq) , which also offers itself as a natural 
choice of dynamical guide? To the extent that our phenomenological scheme 
proved to be robust against such modifications of its ad hoc elements, we 
might feel more confident that it adequately realized the underlying idea of 
A fluctuations of typical magnitude 1/y/V. Investigating this question of 
"structural stability" is the purpose of the present paper. 

2 The "Mixed Equation" and Structural Stability 

The nice results of [3] were obtained from HCEq (2). To what extent would 
they persist if we used instead an arbitrary linear combination of HCEq with 



4 That is, we can retain both Einstein equations if we treat the cosmological term as a 
perfect fluid, but make the altered identifications, pA = A, pa = —A — A/3H. 
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AccEq? 

The first thing to notice in response to this question is that not all linear 
combinations make sense, unfortunately, because AccEq (3) is inherently 
unstable, even with a non- fluctuating A (indeed, even for A = 0). In the 
latter setting, equation (3) is actually redundant since it is (up to a fac- 
tor) just the time derivative of equation (2). Conversely, solving (3) with 
initial data that exactly satisfies (2) is equivalent to solving (2), as is well 
known. But if one's initial data fails slightly to satisfy (2), the failure grows 
with time. This makes (3) unsuitable for numerical solution, even when 
A = 0. Still less would one expect it to be suitable for a stochastic model 
like ours, where random deviations from the constraint submanifold are con- 
stantly being introduced dynamically, even forgetting about the round-off 
and discretization errors of one's differencing scheme. 5 

The instability in question is brought out clearly if one re-expresses 
AccEq in terms of the variable D = 3H 2 — A — p, which measures the 
degree to which the HCEq fails to be satisfied. Doing this yields for the 
time-dependence of D 

- + D = -±, (4) 
3H 3H v 1 

or 

a 3_ (a 3 D) = __ , ( 5 ) 

while HCEq itself is, of course, equivalent to 

D = . (6) 

For a A which doesn't fluctuate, say for A = 0, the right hand side of 
equation (5) vanishes, and the time dependence of D is particularly simple: 

a(r) 3 D(r) = constant, (7) 



or 

DW ^)Vo). (8, 

Noting that by definition 3H 2 = A + p + D, and recalling that an a~ 3 
scaling is precisely that of "dust" (pressureless matter), we can interpret 
(8) as saying that any deviation from D{tq) = acts as a fictitious source 
of dust introduced at time tq. In a radiation-dominated cosmos, such a 
term will ultimately swamp the genuine stress-energy, no matter how small 
it starts, since p ra diation scales as a -4 . At that point, the relative error 
D/p reaches O(l) and the behaviour of the solution becomes very different 
from what it would have been had the error not been introduced. If Dq 

5 The instability of the "constraint submanifold" has often plagued attempts to solve 
the Einstein equations numerically. 
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happened to be negative, the cosmos would re-contract, even though in the 
true solution it would have continued expanding forever. (Notice also that, 
for numerical purposes, the instability is much worse than (8) would make 
it seem. Because a varies by a factor of 10 30 or so over the course of a 
simulation, it is impractical to use r as time-parameter. Rather, something 
like log r must be used, and with respect to such a time, the instability will 
be exponential.) 

For a time- varying A we find instead of (8), 

D <H^) 3 W>'(^) 3 <-AM>, (.) 

which says, in effect, that a variation dA introduces into the effective energy 
density a fictitious dust contribution of — dA (as if an amount dA of A had 
"turned to dust"). Although this expression is less easy to analyze than 
(8), it is hard to believe that the fluctuations in A under the integral sign 
could avoid exciting the instability. Indeed, simulations we performed with 
a fluctuating A exhibited an instability in this case that was as bad as or 
worse than that for A = 0. 

As a check on this conclusion, we also simulated AccEq another way, 
which is perhaps worth reporting here since it brings out the fact that our 
scheme can be interpreted as a modification to the "equation of state of A" , 
rather than as a change to the Einstein equations. In this re-interpretation, 
we write the Einstein equations as 

3H 2 = PA + p, ^ + H 2 = -( PA +p) . 

The scheme based on (2) is then equivalent to the equation of state, 

PA = A , PA = -A- , (10) 

according to which pa depends on both A and its time-derivative, 6 while the 
scheme based on (3) is equivalent to the equation of state 

PA = -A, pa + pa/3H = A, (11) 

according to which pA is not even a local (in r) function of A. As expected 
(since the two approaches are equivalent modulo numerical errors) simula- 
tions with this scheme exhibited the same instability as did direct simulation 
of equation (3). 

Now let us extend this stability analysis to the equation formed as an 
arbitrary linear combination (with coefficients v, p) of HCEq (2) with AccEq 

6 An equation of state like this was considered in [12]. Notice that, strictly speaking, 
dA/dr is infinite because A is (in the continuum limit) a function of Brownian type. 
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(3). Expressing the resulting equation in terms of D yields the corresponding 
linear combination of (6) and (4), namely 

»£j + (n + v)D = 0, (12) 

where we have again taken A = 0. Instead of (7), we find now 

„ constant 

D = —^r~ > (13) 

where b = (// + v)j[i = 1 + v j p. In a radiation-filled cosmos, we therefore 
need 36 > 4, in order that the error remain small. The limiting case is 
thus b = 4/3 or p = 3v. Interestingly, this corresponds exactly to the 
trace equation, TrEq (given by v = 1, p = 3, and equivalent to G% + 4A = 
T£), which accordingly is "marginally unstable" in the structural sense. 8 
The trace equation (TrEq) and the Friedmann equation (HCEq) are then 
at the two extremes of the stable range, and one may expect that every 
other (stable) combination will behave in a way intermediate between these 
two extremes. This makes TrEq all the more interesting, supplementing its 
intrinsic interest in relation to conformal variations of the metric. 

In the rest of this paper, we will refer to the above linear combination 
of AccEq with HCEq as MixedEq, and we will limit p/v to its stable range. 
In terms of conformal time, MixedEq is given by 

a - 3is (a') 2 a 3 a 3 
pa = — — + — [y ptotal - V Ptotal) + y (A* + v ) A (14) 

with 

Ptotal = Pradiation T" Pmatter j 

and 

_ Pradiation 
Ptotal — g • 

Here p ra diation /matter is the energy density in radiation/pressureless matter, 
and a prime represents differentiation with respect to the conformal time, 
n = J dr/a. For u = 3v we obtain the TrEq written as 

a" 4 A + p ma tter M ^ 

6 • (15) 

(Notice that p ra d drops out of TrEq, as it had to given the conformally 
invariant nature of radiation. For a radiation filled universe with Pmatter — 



7 For completeness we mention that if later on the universe becomes "matter dom- 
inated", since in that era mass density dilutes as a -3 , the inequality can be loosened 
slightly. 

8 At the other end of the stable range, HCEq betrays no hint of instability when con- 
sidered on its own. But adding even a small amount of AccEq increases its differential 
order, whence it is not strange that such a change can destabilize it when the sign of the 
addition is unfavorable. 
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A = 0, (15) yields simply da/dr/ = constant.) We will see in the next section 
that equation (14) with < [i < 2>v is indeed stable and produces results 
similar to those obtained in reference [3] from HCEq. 



3 Results of simulations 

Computer simulations of MixedEq (14) in the stable regime, ji < 3is, produce 
results similar to those obtained in [3] from HCEq. Figure 1 displays plots 




Figure 1: Tracking behaviour exhibited by |A| for various values of the 
model-parameters \xjv and a. 

of the energy density versus the scale factor in some of these simulations, 
and the tracking behaviour alluded to earlier is clearly visible. The top two 
diagrams cover the whole period from the Planck time to the present epoch, 
while the bottom two diagrams focus on the transition from a -4 scaling 
(radiation domination, shown by the dotted line) to a -3 scaling (matter 
domination, shown by the dashed line) of the ambient energy density. One 
sees that the effective energy density in A (solid jittery line) switches its 
scaling as well and follows the total effective energy density (wiggly dotted 
line) modulo fluctuations. 

One of the benefits of using the MixedEq (i.e., any of the stable combi- 
nations with non-vanishing jjl) is the ability to accommodate larger values 
of the parameter a. With HCEq almost all of the simulations terminate 
prematurely when a is too big, the reason being that when a is large, large 
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Figure 2: Recollapse is typical of evolution with MixedEq (or TrEq), and 
prevents the cosmos from reaching its present size when a is too large. In the 
example shown a = 1/5. The energy density in matter and radiation (solid 
line p) moves from left to right as the cosmos expands, then retraces its path, 
heading toward the final singularity. In contrast, |A| (dashed line) tracks p 
while the cosmos expands, but remains essentially frozen as it contracts. 



fluctuations in A are bound to happen, including negative ones. If such a 
fluctuation overwhelms p in (2), it renders further evolution impossible by 
making H imaginary. There are intriguing suggestions of how to keep going 
in such a circumstance, some of which were discussed in [3], but a good 
understanding of their status is still lacking. 9 

On the other hand for MixedEq (14) with < p < 3z/ a large negative 
fluctuation in A simply means that a becomes negative. 10 This can even- 
tually turn an expansion into a contraction (and vice versa for a positive 
fluctuation), but it poses no problem of principle and it lets the simulation 
continue. Thus we have been able to simulate as high as a 2 = 1/3 with 
MixedEq, whereas a 2 > 1/50 was not viable in [3]. The use of larger values 
of a has another, important benefit as it naturally results in larger final 

9 The most obvious suggestion is simply to change the sign of a at that point, but it 
seems that this alone cannot cure the problem in general. 

10 Actually a" becomes negative but this implies that a is also negative as can be seen 
from the equation a = — . 
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values of U\. We will discuss this further in connection with figure 6. 

Although MixedEq permits us to evolve through a turning point, thereby 
affording a larger latitude in the choice of a, simulations with a > 1 still do 
not reach the present value of the scale factor (or equivalently the tempera- 
ture) in any significant number before re-collapsing. Thus, the problem that 
affected HCEq [3] shows up in another guise with MixedEq, albeit at larger 
values of a. The biggest a that reached the present epoch was a/1/3 0.58, 
and out of 10 million runs with this a only 25 did so. (The first graph in fig- 
ure 1 shows one of these simulations) . In order to explain what goes wrong, 
recall that the magnitude of A in our model decreases with the past volume 
(which, given our ansatz, is constant on each hypersurface of homogeneity). 
As long as the cosmos is expanding, the absolute magnitude of the A term 
remains comparable to the energy density in radiation (this being the rel- 
evant energy component at early times). 11 But when the cosmos begins 
to contract, the density of radiation, being proportional to a -4 , increases 
sharply. And since the past-volume keeps on accumulating, after a while the 
A fluctuations are too small to reverse the sign of d, and the probability for 
the universe to recover from the collapse practically vanishes. The graph of 
Figure 2 illustrates one such case. 




o 

1 

0.2 - 



0.05 0.1 0.5 1 1.5 2 2.5 2.9B7 2.998 2.999 3 3.001 3.002 
fi/v 

Figure 3: The fraction of runs that reach the present time, plotted as a 
function of fi/v, for a~ 2 = 50. The HCEq corresponds to fj,/v = 0, the 
TrEq to n/u = 3. 

As an indicator of how our model reacts as we move away from the HCEq 
(/j, = 0) into the MixedEq (0 < fx < 3v) and then reach the other extreme at 
the TrEq (fi = 3^), we studied the fraction of times our simulations reached 

n This of course is the tracking behaviour that solves the "Why Now?" puzzle. 
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0.1 0.2 0.30.5 1 1.5 2 2.5 2.97 2.98 2.99 3 3.01 3.02 
fi/v 



Figure 4: Same as figure 3 but for a 2 = 30. 

the present value of the scale-factor. For a 2 = 1/30,1/50, and 1/70, we 
simulated HCEq, TrEq, and MixedEq with values of fi/v as small as 0.002 
(just after the HCEq) and reaching as close as 2.99998 to the TrEq (p,/u = 3). 
Figures 3 and 4 plot this fraction as a function of p,/u for a 2 = 1/50, 
and 1/30, respectively. In both plots, one observes a central "plateau" in 
which most realizations reach today's scale-factor, flanked on each side by 
"escarpments" in which the fraction of such realizations falls precipitously 
as one approaches the unstable range. Close to the TrEq the dependence 
seems to be exponential, so in figure 5 we have expanded that part of the 
curve for a 2 = 1/50 to reveal the continuous transition to p/u = 3. 

Figure 6 shows some histograms of the final values of f^A (called Q^ nal ) 
for HCEq with a 2 = 1/60 and for MixedEq with a 2 = 1/30, and for various 
pjv values. It can be seen that the values of f2 A maZ are nicely peaked at 
zero and fall off on either side as expected. The first thing to notice is 
that S7 A ma ' goes as high as 0.99. About 5 percent of the time, it lies above 
0.5, and about 2 percent above 0.7. The positive values of VL^ nal slightly 
outweigh the negative ones, the probable reason being that a universe which 
spends more of its time with positive A is more likely to reach the present 
epoch without recollapsing. In general the positive values seem to make up 
fifty to sixty percent of the total. One more thing that catches the eye is 
that the tail on the negative side is longer than on the positive side. This 
merely reflects the asymmetry built into the definition, f^A = Pa/(pa + Pm), 
in consequence of which J7a normally 12 is bounded by +1 on the positive 

12 It can happen that pa + p m becomes negative such that J7a > 1- However the re- 
contraction that accompanies a negative net effective energy density enhances the matter- 
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2.9992 2.9996 3 



Figure 5: The end part of the curve in figure 3 magnified about twenty 
times. 

side, while ranging to — oo on the negative side. (A few percent of the JIa 
fell below —1, but we did not plot them on the histograms.) 

4 Conclusions 

This study is a follow-up to reference [3], where we recalled how the causal 
set idea leads to fluctuations in the cosmological term and showed that if 
the mean value of A is taken to be zero then these fluctuations have the 
potential to account for the presently observed Oa- We also raised some 
concerns about the model, namely: 

• The nice results of [3] were obtained from the choice of equation (2) 
as "dynamical guide", and it was uncertain whether (3) or a linear 
combination of the two equations would behave similarly. 

• It was not clear how to continue the evolution when the total effective 
energy density became negative. 

• The parameter a that governed the magnitude of the fluctuations had 
to be five to ten times smaller than unity, and this looked like a mild 
fine tuning problem. 

• Rarely did the model produce a final value for Q\ large enough to 
accord with present observations. 

density, and Qa rapidly returns to its "normal" range, either for this reason or because A 
fluctuates back to a positive value. 
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a* = 60, |Vv = zero (HCEq) 



= 30, n/v = 2.7 




1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 

linal 



Figure 6: Histograms of the final values of for fi/u = 1.5, 2.7 and 3 
(a = ^/l/30). Also shown for comparison is a histogram for fi = and 



The present work throws light on all of these issues. We take them in turn. 

• It turns out that the second order equation (3), referred to in this 
paper as the AccEq, is inherently unstable and cannot replace HCEq 
(2) in the model of [3]. However a large range of convex combinations 
of the two equations is stable, specifically combinations of the form 
v x HCEq + fi X AccEq with < fi < 3u. In our simulations, all 
of these stable combinations exhibited the same tracking behaviour 
as was obtained in [3] from the purely first-order HCEq (p = 0). 
Our question about "structural stability" is therefore answered in the 
affirmative. 

• When the net effective energy-density of a universe (evolving under 
one of the stable mixtures of HCEq and AccEq) becomes sufficiently 
negative, an expanding universe starts to contract. In principle, it can 
re-expand if the net effective energy-density becomes positive again, 
but this happens only rarely. In any case, there is for fi > 0, no 
problem in following an expansion and re-contraction all the way down 
to a singularity of infinite density. 

• Whereas in [3] with the HCEq, we were only able to simulate as high as 
a 2 1/50, we have been able with MixedEq to go as high as a 2 = 1/3, 
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albeit 1/30 is a more realistic upper limit if the universe is to reach its 
present size with an appreciable probability. 

• Since with MixedEq we can handle larger values of a (with their con- 
comitant larger fluctuations) , we can end up with larger final values of 
Qa . For a ~ 1/30 we can account for the present observational value 
of Qa in about 2 percent of the simulations. 

It is important to realize, however, that our variable cannot be compared 
directly with observational parameters quoted in the literature which assume 
that A is constant (or has an "equation of state" with constant w = pa/ PA )• 
When instead, A is fluctuating, the comparison with observation must be 
derived anew. To that end, one might for example use our simulations to 
construct fictitious data sets of supernova luminosity and redshift, one for 
each run of the simulations. One would then feed these data sets into one 
of the algorithms people have used to obtain the quoted values of A and w, 
and one would ask how often the resulting values came near to the quoted 
ones. 

Can the accelerating Hubble expansion be traced to a fluctuating A and 
can such fluctuations be understood as a nonlocal and quantal residue of an 
underlying spatio-temporal discreteness? 

This long-standing idea has still to be embodied in a fully fledged phe- 
nomenological model, but an encouraging start was made in [3]. The model 
of [3] involved some ad hoc choices, however, together with an artificial 
restriction to spatial homogeneity. In this paper we have not dealt with 
the latter shortcoming 13 but have concentrated on the ambiguity inherent 
in the manner in which [3] implemented the idea of a varying A. So far, 
the evidence is that none of the qualitative features of the model depend 
on how that ambiguity is resolved, and overall our results sustain the pic- 
ture developed in [3], according to which the fluctuating A tends to remain, 
throughout the phase of cosmic expansion, in rough equilibrium with the 
ambient matter density. 

There remains, though, a kind of tension between the magnitude of the 
fluctuations (as reflected in the parameter a) and the continuation of the 
expansion. If a is too big, the negative fluctuations will tend to terminate 
the expansion; if it is too small, the positive fluctuations will be unable to 
account for the current value of A. If one chooses a appropriately, the two 
competing effects can be balanced 14 fairly well, but some discrepancy always 
seems to remain. Indeed, it is likely a generic feature of the models based 
on MixedEq that for any choice of a, the cosmos given long enough, will 
eventually recollapse due to a negative fluctuation. 

13 The consequent risk to the model was pointed out in [3] and emphasized further in 
[18]. 

14 Order of magnitude estimates suggest that this balancing is possible only in 3 + 1 
dimensions, as pointed out also in [18]. 
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If this is so, we hope we may be excused for speculating further that 
cycles of expansion and contraction would succeed each other indefinitely, 
their characteristic lifetime depending on the value of a as it emerges at 
the start of each new cycle. In this way, the "Tolman-Boltzmann" scenario 
of [19] would have acquired a possible basis in quantum gravity. Although 
that scenario emerged from a non-quantum dynamics for causal sets, namely 
that of classical sequential growth, it illustrated how the most striking fea- 
tures of our universe — its approximate spatial homogeneity and isotropy 
- might have emerged dynamically over the course of repeated expansions 
and collapses. There the collapses occurred via rare statistical fluctuations 
following exponentially long periods of stasis. Here, the collapses would be 
more dynamical in character, initiated by quantal fluctuations in A, or in 
other words fluctuations in the form of the gravitational field equations. 

This research was supported in part by NSERC through grant RGPIN- 
418709-2012. Research at Perimeter Institute for Theoretical Physics is 
supported in part by the Government of Canada through NSERC and by 
the Province of Ontario through MRI. 
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